Relativistic effects in quantum walks: Klein's paradox and Zitterbewegung 
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Quantum walks are not only algorithmic tools for quantum computation but also not trivial models 
which describe various physical processes. The paper compares one-dimensional version of the free 
particle Dirac equation with discrete time quantum walk (DTQW). We show that the discretized 
Dirac equation when compared with DTQW results in interesting relations. It is also shown that two 
relativistic effects associated with the Dirac equation, namely Zitterbewegung (quivering motion) 
and Klein's paradox, are present in DTQW, which can be implemented within non-relativistic 
quantum mechanics. 
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Since the_papers of Aharonov et al. and Farhi and 
Gutmann quantum walks have been deeply inves- 
tigated in hope to find faster algorithms (Refs. 0, ^ 
and references therein). Despite their contribution to 
the quantum information processing, quantum walks are 
themselves very interesting physical systems worth be- 
ing studied due to effects from various fields of physics 
like quantum chaos 0| and solid state physics 0. Re- 
cently, it has been shown that discrete time quantum 
walk (DTQW) resemble the one-dimensional free particle 
Dirac equation dli. Actually, one has to keep in mind 
that the idea of quantum walks goes back to Feynmann 
and Hibbs 01 who considered a discrete version of the 
one-dimensional Dirac equation propagator. In this pa- 
per we study differences and similarities of the two mod- 
els and show that two effects associated with the Dirac 
equation, Klein's paradox and Zitterbewegung, occur in 
DTQW. This is an important result leading to the fact 
that the Dirac equation might be simulated via DTQW. 

The one-dimensional free particle Dirac equation has 
the form 
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where m is a mass of the particle, and at and aj have 
to obey af = — I and ataj + Ujai = 0, as we want 
Eq. to be Lorentz invariant. We assume c = % — 1 
here and throughout the paper. The most natural way is 
to take Ui and Uj to be the two distinct Pauli matrices. In 
this case the wave function has two components ^{x, t) — 
{4>u{x,t),ipd{x,t)) . The eigenvalues of Eq. ^ a.re E± = 



where k is the eigenvalue of the momentum 



operator. The striking feature of the Dirac equation is 
that it gives positive and negative energy solutions and 
that there is a gap of forbidden energy region —m < E < 
m. 

Now, let us define DTQW. The step operator of one- 
dimensional DTQW is given hy U = TC, where 



(2) 



is the conditional shift operator indroduced by Aharonov 
et al. 0, and C is a two-dimensional coin operator we 



assume in the form 
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In Eq. Px, <^z and I are the x component of momen- 
tum operator, the z component of spin vector operator 
and the size of one step respectively. In Eq. 6* is the 



angle of rotation about n 



I Tirf , Til, , ^7. 



which is an 



arbitrary unit vector, and a — {o'x,o'y,az)'^ is the spin 
vector operator consisting of three Pauli matrices, there- 
fore a± = cos I ± zn^ sin I and /3± = {irix ± rty)sin|. 
The state of the walk is described by the two component 
wave function ip{x,t) = {tjj+{x,t),tp-{x,t))'^ , like in the 
one-dimensional Dirac equation, and from Eqs. (01 and 
one can obtain a recursive formula 

V'i {x, t + At)^ a±il;±{x Tl,t)+ l3±ij+ {x Tl,t), (4) 

where At is the time of one step. Taking n — 
(^,-^,0)-^ and 9 = TT one gets the well studied 
Hadamard walk 

To calculate eigenvalues of U we have to realize that 
it commutes with the translation operator and therefore 
with the momentum operator. Thus, the eigenvectors 
are of the form |A) = N{^^ e''''^\x)) (8) |c), where |c) is 
a position independent state of the coin and A'^ is a nor- 
malization constant. Together with the coin operator C, 
Eq. ||2Jl, this gives eigenvalues of the form 



'^it (fc) — cos - cos k + Uz sin - sin k ± 



i Y 1 — (cos - cos k + Hz sin - sin fc)^. (5) 

It is easily verified that A5-A± = 1 and one can write them 
as A±(fc) ^ e"*'^±('=) = e-'-^±('=)'^*, therefore we may as- 
sume that A_|_(A:) and A_(fc) correspond to positive and 
negative energies respectively. On the complex plane, 
X±{k) are confined to two symmetrically inverted cones, 
one corresponding to positive while another to negative 
energies (see Fignj. The dilation angle of the cones 
and the angle between the real axis and the line passing 
through the middle of the both cones depend on 6 and . 
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FIG. 1: Eigenvalues of DTQW evolution operator on the 
complex plane. +E and —E denotes the regions of positive 
and negative energy respectively. 

For the Hadamard walk, the both angles are 7r/2. The 
region of complex unit circle not covered by A-i-(fc) might 
be thought of as a forbidden region. One sees that the 
above properties make DTQW and the Dirac Equation 
alike. 

The continuous limit of DTQW, both in time and 
space, has been studied in Refs. It was shown 

that the continuous version of DTQW gives the evolution 
similar to the one described by the one-dimensional Dirac 
equation and that both models give the typical two horn 
probability distribution for initially highly localized wave 
packets. Aslangul 10] has studied the model with a coin 
of an arbitrary dimension with the coin operator being 
rotation of spin about y axis, so the exact comparison 
with the Dirac equation cannot be made, although the 
models seems alike. It is also important to notice that 
the coin operator makes DTQW to be more general than 
the Dirac equation. In our case it can be an arbitrary 
2x2 unitary matrix. Due to this fact, the continuous 
version of DTQW must not always be Lorentz invariant. 
This is the case of the Hadamard walk. 

Now, we will do something opposite to what was done 
in Refs. "t", 'l^. Let us discretize Eq. and see if it 
gives Eq. Q . We take ct^ and aj to be the z and y Pauli 
matrices respectively and use the same time symmetric 
procedure as in Ref . Tl'l . We obtain a recursive formula 
that conserves probability to the second order in At: 

^ {x, t + At) = ^^^V^ {x T Ax, t) T i (2^' 

^ ^ (6) 

where 7 — mAi/2. This looks similar to Eq. but 
to make it more alike one should shift the right hand 
side of Eq. © by / = Ax/ 2 with the conditional trans- 
lation operator T, Eq. It means that DTQW is the 
discretized Dirac evolution viewed in, somehow strange, 
conditional moving reference frame. Note, that this shift 
does not occur when one considers the continuous limit 
of DTQW. This is due to the fact that in the continu- 
ous limit the differences between U = TC and U = CT 
vanishes because C and T act simultaneously. In the dis- 



FIG. 2: Klein's paradox. The forbidden region —m < E < m 
is shifted by a potential V. The particle with positive energy 
E is transferred to the negative energy region. The wave 
function of particle is not damped in the new region. 

Crete regime it is important wether we act with TC or 
CT. One can also easily check, that in the case U — CT, 
the left hand side of Eq. ^ should be shifted with the 
operator in order to resemble the DTQW recursive 
formula. Also, choosing ai and aj to be different than 
here, gives the same effect. One also needs to remem- 
ber that since DTQW is by definition discrete in space 
and time, it should be rather compared with the discrete 
Dirac equation than with its continuous more common 
version. It is due to the effects like Bloch oscillations 
caused by discrete space and perhaps due to the other 
effects, yet unknown, related to discrete time. 

Klein's paradox and Zitterbewegung are effects that 
were first observed for the Dirac equation (for reference, 
look for example Ref. ^3)- The first one is transfer of 
the particle with energy E from the region with zero po- 
tential to the region with potential V > V — m > E 
(see FigEJ. In relativistic case the wave function is not 
damped inside the potential region, unlike the solutions 
of the Schrodinger equation. The second effect is the 
rapid oscillation of a wave packet due to interference of 
its positive and negative energy components. 

To observe Klein's paradox in DTQW, one has to 
introduce a potential V . For simplicity, we take the 
step potential: all over the region x > a the potential 
is Vq and in the region x < a the potential is zero. 
Adding the uniform potential simply shifts the energy 
E' = E + Vo and in the exponential form it may be writ- 
ten as e"'^'^* = e-'(^+^o)^* = e'^^^^e^^. From now 
on, we will identify the potential with e*'''. This will be 
very helpful to describe the existence of Klein's paradox 
in DTQW in geometric way. 

The eigenvalues without the potential are presented in 
FigEla. In the presence of the potential, the eigenvalues 
are rotated by the angle ^p, see Fig|31 b. Imagine that 
DTQW starts in the region without the potential and 
that there is a wave packet heading for the potential step 
giip_ The transmitted part of the wave packet may be- 
have in two distinct ways: either move further without 
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FIG. 3: Eigenvalues of DTQW evolution operator on the 
complex plane: a) eigenvalues without potential, b) the same 
eigenvalues in the presence of a potential e**^, c) intersection 
of the eigenvalues of two regions, one without potential and 
one with the potential e"^, giving non decaying wave func- 
tions (dark grey), d) intersection of the eigenvalues of two 
regions, one without potential and one with the potential 
g2(¥)+7r)^ giving non decaying wave functions and leading to 
Klein's paradox (dark grey). 



damping, or start to decay exponentially fast. Whether 
it is the first or the second case, depends directly on 
the eigenfunction decomposition of the wave packet. Let 
{A}o denote the set of the eigenvalues of DTQW with no 
potential and {A},^ be the set of the eigenvalues in the 
presence of the potential. Moreover, let {X}wp C {X}o be 
the set of the eigenvalues which correspond to the eigen- 
functions that contribute to the initial wave packet. By 
calculating {Aj^pnlA};^ one obtains the part of the wave 
packet that will not undergo damping. Now, let us study 
how this depends on tp. For small ip the wave packet 
is not damped if its energy is higher than the potential 
step, similarly as in non-relativistic quantum mechanics, 
see Figl^l c. As (p rises, more eigenvalues fall into the 
forbidden region and the corresponding eigenfunctions 
are damped, until if is large enough and the eigenval- 
ues from the cone {E- ) intersect with the eigenvalues 
of the rotated cone E-{E^) (FigOld). This is exactly 
Klein's paradox. We should also point out that for spe- 
cific coins, namely for the coins which make the eigen- 
values cone dilation angle less or equal to 7r/2, one can 
choose a potentials that damp all eigenfunctions. The 
Hadamard walk is always damped in the potential region 
if </? = 7r/2. This might be useful if one would like to con- 
fine DTQW to some region. For example, it is possible 



to study the quantum walk on the line with one or two 
reflecting walls. 

Zitterbewegung oscillation are expected to occur in 
both position and velocity. To measure this effect, we 
calculate how the position changes with the one step 
AX ^ X{t + At) - X{t). This is the discrete version 
of velocity. In the Heisenberg picture this might be writ- 
ten as AX = WXU ~ X, where X = T,x^\^)i^\ and 
U — TC, as before. This leads us to the position inde- 
pendent coin operator 
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To calculate its time dependence, we perform the discrete 
Fourier transform on the initial state and use the momen- 
tum representation to derive the expectation value of AX 
as a function of time. Calculation of DTQW probabil- 
ity distribution via Fourier transform was presented by 
Nayak and Vishwanath . The general form of the time 
dependent wave function in momentum representation is 

^(fc,i) = ^/,(fc)(A,(fc))*|c,(fc)), (8) 

where f±{k) are any square integrable functions obey- 

ing J^:dk{\U{kr + \f-{kr) = 1, and \c±{k)), 
(czfz{k)\c±{k)) = are the coin states corresponding 
to the eigenvalues A±(fc). Since {c-{k)\AX\c+{k)) = 
{c+{k)\AX\c-{k))* — g{k), the expectation value of 
AX{t) may be written a,s A + B{t), where A is a constant 
corresponding to the uniform velocity of the wave packet, 
and B{t) is a time dependent term 



B{t) 



dk (/;(fc)/_(fc)g(fc)(A;(fc)A_(fc))* + c.c.) 



(9) 

The above vanish if the initial wave packet consists only 
of the positive or the negative energy eigenfunctions. As 
before, one may write the eigenvalues in the exponential 
form A±(fc) = e-*'^±('=) = g-^-^i^^)^*, so the time depen- 
dent part under integral, Eq. Q, yields 



(A;(fc)A_(fc))* -e'^^+W-^-W)*, 



(10) 



which has the form of oscillations. If we take the initial 
wave packet largely spread over the position space (for 
example broad Gaussian packet) and multiply it by the 
factor e*'^°^, the corresponding momentum wave packet 
would be mainly localized around fco- In this case Eq. Q 
can be approximated by 



Bit) w Bcos{{E+{kQ) - E^ika))t) 



(11) 



where _B is a time independent amplitude of oscillations. 
Of course, time is discrete and t = nAt. 

Zitterbewegung is very well visible if one chooses the 
standing wave packet, what in the case of DTQW does 
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not always mean fco = 0. It is visible if one calculates the 
group velocity 

Without loosing generality we assume At = 1 and obtain 
n, sin ^ cos k — cos ^ sin k 



v±{k) = ± 



1 — (cos I COS k + riz sin | sin fc)^ 



(13) 



Since the wave packet represents particle which does not 
move, {E+{ko) — E-{kQ))At has the meaning of the mass. 
This result is closely related to the one of the Dirac equa- 
tion with the oscillation frequency approximating 2m. 
For the Hadamard walk, v±{ko) — for fco = tt/2 and 
A±(fco) = thus (E+iko) - E^{ko))At = n/2. 

The presence of Klein's paradox and Zitterbewegung in 
DTQW leads to the bunch of questions. First of all, are 
the two effects truly relativistic, since DTQW might be 
implemented with non-relativistic quantum mechanics? 
The implementations of DTQW has been presented for 
vario us phys ical systems, from optical latices to trapped 
ions 111 [IE 111 111 El 111, just to point some of them. 
One may even think of the Stern-Gerlach experiment as 
of a quantum quincunx with the magnetic field along 
X and z axes at even and odd time steps respectively. 
The implementation schemes might be different, but it 
is worth to notice that none of them has anything to do 
with the Dirac equation. Of course one may wonder if it 
is right to consider such concepts like spin and talk about 
non-relativistic quantum physics, but the fact is that we 
may include spin in the non-relativistic Schrodinger equa- 
tion without bothering about its origin. Moreover, this 
equation would correctly predict the behavior of a real 
system. 

In conclusion, effects associated with relativistic quan- 
tum mechanics were shown to appear in non-relativistic 
models. We derived formula for Zitterbewegung oscilla- 
tions, Eq. Ulljl . and presented in graphical way the nature 
of Klein's paradox. DTQW was also compared with the 



discrete Dirac equation and it was shown that the two 
models are related to each other and that one may go 
from one to another by changing the reference frame. 
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